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ABSTRACT 

We discuss one-loop electroweak corrections to the decay width of toponium. We calculate 
the energy-level shifts by expanding around the solution of the Bethe-Salpeter equation in the 
instantaneous approximation, in analogy with the positronium case. We show that first-order 
electroweak effects are suppressed by at least four powers of the strong coupling constant, and 
are therefore negligible compared with QCD corrections. The calculation is manifestly gauge 
invariant and takes into account the contributions to the decay rate due to both Coulomb 
enhancement and phase space reduction effects. 
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The threshold production of a tt pair at future e^e~ coUiders has attracted a lot of atten- 
tion because the measurement of its cross-section will permit a simultaneous determination 
of the mass of the top quark nit, its decay width Tt, and the strong coupling constant as{mt) 0. 
The threshold regime is dominated by tt resonances which enhance the cross-section. These 
bound states differ from those of the other heavy quarks (b and c) , because rrit is so large that 
the dominant decay mode of the top quark is the weak decay t — > W^b. Thus, the decay width 
of a top quark is 

r. - ^< . (1) 

For rrit ~ 174 GeV , we obtain a decay rate of the order of 1 GeV, which exceeds Aqcd, 
and permits the use of perturbation theory, because the tt pair will decay before it has time to 
hadronize. Moreover, the quarks are produced with little kinetic energy, so that non-relativistic 
techniques are applicable. Thus, the toponium energy levels are well approximated by the 
Bohr levels of a QCD Coulomb-like effective potential. To calculate the cross-section, one can 
determine the appropriate Green function which obeys a non-relativistic Lippmann-Schwinger 
equation with the QCD effective potential p|. 

An important parameter entering the Lippmann-Schwinger equation is the decay width of 
toponium, Fjf. In general, it can be a function of the momentum of the top quark, and its 
precise form is crucial for an accurate theoretical prediction of the cross-section cr(e+e~ — > tt). 
To zeroth order, 

r^f = 2r, . (2) 

This relation is modified by three effects: (a) time dilatation (the top quark lives longer in the 
center-of-mass frame), (b) phase space reduction (due to the binding energy of the quarks), and 
(c) Coulomb enhancement (the b and b quarks should be described by Coulomb wave functions 
rather than plane waves) The first two effects can be easily accounted for and they lead 
to a significant modification of Eq. (§). The third effect is much harder to take into account. 
As was shown by Kummer and Modritsch, it cancels the second effect to 0(0^), thus leaving 
time dilatation as the only effect modifying Eq. (|]), in analogy with the case of muonium in 
nuclei [Q. We therefore obtain 

Tu = 2Tt{l-{pVm^)f' , (3) 
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which is only a small modification of Eq. (|^). To arrive at this result, Modritsch and Kum- 
mer |^ calculated one-loop graphs with insertions of bound-state wave functions, as prescribed 
by the rigorous Bethe-Salpeter formalism ||^. The cancellation of the two effects (Coulomb 
enhancement and phase space reduction) was due to non-trivial cancellations between Green 
functions, which went well beyond the cancellation of gauge-dependent pieces implied by Ward 
identities. 

Our purpose here is to extend the results of ref. to higher orders in the strong coupling 
constant. We shall show that electroweak corrections enter at an order higher than o (a^), and 
therefore we need second-order bound state perturbation theory in order to calculate them. 
They are due to corrections to the vertex function, which do not contain a contribution to 
the magnetic moment of the top quark. This is in contrast to the positronium case, where 
magnetic moment effects are o (a^). The discrepancy arises from the fact that W only couples 
to a left-handed current. 




Figure 1: The homogeneous Bethe-Salpeter equation. 



We start with a brief review of the Bethe-Salpeter formalism in order to fix the notation. 
We wish to calculate the energy levels of tt bound states. They are poles in the four-point 
amplitude describing tt scattering. We shall only consider scattering in the t-channel, because 
the annihilation diagrams are suppressed 0]. We shall also neglect photon and exchanges 
because they are small effects compared to a gluon exchange. A Higgs exchange will also be 
neglected, but depending on the mass of the Higgs boson, it can have an appreciable effect. 
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Figure 2: The zeroth-order contribution to the QCD potential. 
A bound-state wavefunction Xpip) satisfies the homogeneous Bethe-Salpeter equation 

n(^)(p+)n(2)(p_)xp(p) + / ^vM-^ p)xp{p') = o , (4) 

where 11 (p) is the complete inverse fermion propagator, 

n(p)=/-M,-S(p), (5) 

and we have defined momenta 

p 

P± = ^±P- (6) 

Eq. is represented graphically in Fig. [l|. V{p,p'; P) is a potential function which consists of 
the two-fermion irreducible graphs. For our purposes, the mass is complex, 

Mt = mt + iVt . (7) 

To lowest order in and neglecting electroweak interactions, the potential is (Fig. ^ 

Voip, p'; P) = Cf Anas^l^^G^^ip - p'H'^ , (8) 

where Cp = 4/3 is the Casimir operator, and G^^{k) is the lowest-order gluon propagator. In 
the Feynman gauge (omitting group theory factors), 

and the potential Vo{p,p'\P) is independent of P. At threshold, the quarks move with 
non-relativistic velocities and the Bethe-Salpeter equation can be approximated by the non- 
relativistic Schrodinger equation in momentum space, and then solved. To this end, we shall 
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work in the total rest frame in which the overall momentum is P'^ = {E, 0). In the instantaneous 
approximation, the potential becomes 

Vr\p,p'; P) = Cp47cad' \^ ^ 7^'^ ■ (10) 

[p-p r 

If we integrate over p^, we can write the Bethe-Salpeter equation in terms of the wavefunction 

m = I ^Hp) as 

(i/W + H(^^ - E)m = (A? A? - A«aL^)) CAna^J ^.J^^.^ip) , (H) 

where H is the Dirac Hamiltonian and A+ (A_) is the projection operator onto positive (neg- 
ative) energy states. In the non-relativistic limit, this reduces to the Schrodinger equation in 
momentum space 



+ 2M, - e) m = 4.„. / . (12) 

Thus, we obtain the energy levels 



E^ = 2M, - + o{at) , (13) 

which are the Bohr levels of the Coulomb- like QCD potential (plQl). Therefore, the first-order 
QCD correction to the decay rate of toponium is 

in agreement with Eq. (|^) (see 0]). The spherically symmetric S = states are given by 

t,.(p-) ^ iM.C,a.r^'^Jf^ . , ^ . (15) 

where is a polynomial of order n — 1 related to the Laguerre polynomials. For n = 1, we 
have Ci = 16v^27r. 

Higher-order corrections can be systematically introduced by perturbing around the solution 
to the Schrodinger equation (0). The potential to be treated perturbatively is V^ — Vq"'^*. There 
is also a contribution from the disconnected diagrams which are due to the self-energy terms in 
the fermion propagators (Eq. (^), but they can be absorbed in the potential if we make use of 




Figure 3: The first-order electroweak correction to the three-point vertex function contributing 
to the tt decay rate. 

the Schrodinger equation. Thus, according to the Bethe-Salpeter formahsm the first-order 
energy level shift is 

AEn = ($„| Dp{p) (//(I) + _ E^) {V - Vr') {H^''> + H^'^ - E„) Dp{p) , (16) 

where the inner product involves an integral over the four-momentum. H is the Dirac Hamil- 
tonian, and Dp is the product of two free fermion propagators (c/. Eq. (^), which can be 
expressed in terms of the projection operators A± as 

^(1)^(2) 

De{p) = E ^^^2 + pO± - ze)] [E/2 - / ± {Ep - ze)] ' ^^^^ 



where Ep = yp + M/ is the energy of the quark on the mass shell. 

To lowest order, the potential is Vq — Vq"'^* (Eqs. (H) and (p!oD ). This is analogous to the 
positronium case, and produces an o(q;^) shift in the energy levels. The first-order electroweak 
correction is 

V,{p,p';P) = A7rCpasaw{Ai\p+,p'+)G^''{p-pW+ll-'>G^''{p-p')A'^\p.,p'_)^ , (18) 

where p± = P/2 ± p, p'_^ = P/2 ±p', and we have made explicit the electroweak coupling 
constant aw ~ CpM^, where Gp is the Fermi constant and Mw is the mass of the W boson. 
The vertex function A^(p, p') consists of the diagrams shown in fig. |^. It is guaranteed to give 
a gauge invariant contribution by the Ward identity satisfied by the one-particle irreducible 
function, 

{p-prr,ip,p') = uip)~uip') . (19) 



Since we are only interested in first-order corrections, we may replace Mt by its real part 
The contribution of Vi to the energy level shift (Eq. (|16])) can then be written as 

~ 16mt J (27r)4 (27r)4 (p - p')2 + ie {1 + n'^y)^ (1 + nV)" 

where y = Ap'^ /mfCpal and y' = Ap''^ /mfCpal. A simple scaling argument shows that the 
lowest-order contribution to the integral comes from the small three-momentum region. Momen- 
tum insertions contribute additional powers of a^. At low momentum transfer, the three-point 
vertex A^ may be written in general as 

A^ip,p') = Pj^^ik^) + a,,k-:F2{k^) , (21) 

where k = p — p', and a^^ = |[7^ , 7^]. The form factors J-'i and JF2 are regular as /c^ — > 0. 
In the positronium case, JF2 gives an o(a^) contribution to the energy level shift, and is due to 
the magnetic moment interaction. In our case, we need to multiply the gamma matrices by the 
projection operator |(1 — 75), due to parity violation of weak interactions. A straightforward 
explicit calculation shows that the form factor J^2{k'^) vanishes to lowest order in k'^. It follows 
that the three-point vertex is proportional to {p — p'Y (recall that p^ ~p\ = P- ~p'- = P~p')- 

Having established the leading-order behavior of A^, we can now estimate the integral in 
Eq. (pop. As we just showed, A^ contributes a factor {p — p')"^. This factor cancels the gluon 
propagator. Then the integral over pq and p'q can be easily done, because of the respective poles 
in the operators D{p) and D{p'). The resulting expression contains six three-momentum factors 
implying that the integral is 0(0;^). Therefore, the electroweak correction to the decay width is 
negligible. Of course, no conclusion can be drawn regarding the exact value of the electroweak 
correction, because such a high order is beyond the scope of first-order perturbation theory. 

In conclusion, we have calculated electroweak corrections to the decay rate of toponium 
in the dominant mode t —>■ W^b. Our results extend those of Kummer and Modritsch 0], 
who used the Coulomb gauge to obtain a perturbative expansion of the Bethe-Salpeter equa- 
tion. We employed a covariant gauge and perturbed around the instantaneous solution to the 
Bethe-Salpeter equation, in analogy with the abelian case of positronium. We found that the 
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electroweak corrections are suppressed by at least six powers of as- This takes us beyond the 
realm of first-order perturbation theory, therefore we cannot calculate the precise value of the 
correction. It is important to extend our results by including other electroweak effects, such as a 
Higgs boson exchange. The development of an efficient systematic method for such calculations 
involving bound states would be of great interest, in view of the significance of the cross-section 
for threshold tt production at future colliders 
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